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We analyze the low-energy spectrum of a two-electron double quantum dot under a potential bias 
in the presence of an external magnetic field. We focus on the regime of spin blockade, taking into 
account the spin orbit interaction and hyperfine coupling of electron and nuclear spins. Starting 
from a model for two interacting electrons in a double dot, we derive a perturbative, effective two- 
level Hamiltonian in the vicinity of an avoided crossing between singlet and triplet levels, which 
are coupled by the spin-orbit and hyperfine interactions. We evaluate the level splitting at the 
anticrossing, and show that it depends on a variety of parameters including the spin orbit coupling 
strength, the orientation of the external magnetic field relative to an internal spin-orbit axis, the 
potential detuning of the dots, and the difference between hyperfine fields in the two dots. We provide 
a formula for the splitting in terms of the spin orbit length, the hyperfine fields in the two dots, 
and the double dot parameters such as tunnel coupling and Coulomb energy. This formula should 
prove useful for extracting spin orbit parameters from transport or charge sensing experiments in 
such systems. We identify a parameter regime where the spin orbit and hyperfine terms can become 
of comparable strength, and discuss how this regime might be reached. 
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I. INTRODUCTION 

Electron spins in gated quantum dots have been ex- 
tensively studied for their possible use in quantum in- 
formation processing [H-lll- In. this context the main in- 
terest lies in the study of coherent quantum evolution 
of electron spins in a network of coupled quantum dots 
in the presence of external magnetic fields. A double 
quantum dot (DQD) populated by two electrons is the 
smallest such network in which all of the steps necessary 
for quantum computation can be demonstrated. In ad- 
dition, a DQD can host encoded two-spin qubits which 
require less resources for control than the single-electron 
spins in quantum dots. In DQDs, the spins can be ma- 
nipulated exclusively by electric fields in the presence of 
a constant magnetic field, taking advantage of the spin 
orbit and/or nuclear hyperfine interactions [2|. I3-IT3|. 

High precision requirements for the control of spin 
qubits have prompted the detailed study of the inter- 
actions of electron spins in quantum dots. The DQDs 
give experimental access to the coherent spin dynamics. 
Studies of transport through DQDs in the spin block- 
ade regime have been particularly useful for probing 
the electron and nuclear spin dynamics. In this regime, 
charge transfer between the two dots of a DQD can take 
place only when the electrons form a singlet state with 
total spin zero. This allows weak spin non-conserving 
interactions to be studied via char ge sensing [isj or by 
charge transport measurements p^ll6|. even in the pres- 
ence of much stronger spin-conserving interactions. The 
most important spin non-conserving interactions are the 
spin orbit interaction and the hyperfine interaction be- 
tween the electron spins and a collection of nuclear spins 
inside the DQD 



In this work, we investigate the hyperfine and spin- 
orbit mediated coupling between electronic singlet and 
triplet spin states of a DQD in the spin-blockade regime. 
We show that the spin-orbit and hyperfine contributions 
to this splitting can be tuned by a number of parameters. 
We derive an explicit formula that gives this splitting as 
a function of a homogeneous external magnetic field and 
the detuning between the ground-state energies of the 
dots in a DQD. These parameters can be varied in an ex- 
periment. In addition, the splitting depends on the spin- 
orbit coupling interaction and the inhomogeneous nuclear 
Overhauser field, as well as on the dot parameters such 
as the hopping amplitude between the dots in a DQD, 
Coulomb repulsion between the electrons, and the direct 
exchange interaction. Further, we describe how the de- 
pendence of the singlet-triplet splitting on these param- 
eters might be used to extract the intrinsic strengths of 
the spin orbit and hyperfine couplings from charge sens- 
ing measurements in which the DQD is swept through a 
singlet-triplet level crossing in the presence of spin-orbit 
interaction and a fiuctuating nuclear field. 

Recently, it was predicted that the angular momentum 
transferred between electron and nuclear spins in both dc 
transpo rt llSl and Landau-Zener-type gate sweep exper- 
iments 19|, |20| can show extreme sensitivity to the ratio 
of spin-orbit and hyperfine couplings. Our result gives 
an explicit dependence of this ratio on the detuning and 
external magnetic field, thus showing how all regimes can 
potentially be reached. 

The paper is organized as follows. In Sec. |TT]we intro- 
duce a model of a DQD, and describe its energy levels as 
a function of detuning. In Sec. lIIII we find the matrix ele- 
ments of the spin orbit interaction in the space of relevant 
low-energy states. In Sec. IIVI we study the orbital and 
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spin structure of the singlet and triplet states which nom- 
inally intersect for particular combinations of the DQD 
potential detuning and external magnetic field. In Sec.lVl 
we define an effective Hamiltonian which describes the ac- 
tion of the spin-orbit and hyperfine couplings in the corre- 
sponding two- level subspace. Then, in Sec. IVII we study 
the dependence of the resulting singlct-triplct splitting 
on external parameters and show how this dependence 
can be used to extract the spin orbit interaction strength 
and the size of Overhauser field fluctuations from charge 
sensing measurements. In Sec. IVIII we discuss how the 
DQD can be tuned between the regimes of spin-orbit- 
dominated splitting and the hyperfine-dominated split- 
ting. Finally, we summarize our results in Sec. IVIIII 



II. MODEL HAMILTONIAN FOR DOUBLE 
QUANTUM DOTS 

In a DQD, electrons are confined near the minima of a 
double- well potential Vdqd, created by electrical gating 
of a two-dimensional electron gas (2DEG) in, e.g. GaAs, 
see Fig. [T] For the case of a deep potential, we treat the 
two local minima of the double- well as isolated harmonic 
wells with ground state wave functions ipi^2- In order to 
define an orthonormal basis of single particle states for 
building-up the two-electron states of the DQD, we form 
the Wannier orbitals $l and ^ji, centered in the left and 
right dots, respectively (2l| : 

$L,R = , ^ = {^1.2 - gV2,i) , (1) 

where s = {ipi\Lp2) = exp[— (a/ae)^] is the overlap of the 
harmonic oscillator ground state wave functions of the 
two wells, a-Q ~ ^ hjmujQ is the Bohr radius of a single 
quantum dot, hw^ is the single-particle level spacing, and 
2a = Hs the interdot distance. The mixing factor of the 
Wannier states is g = (1 — Vl — s^)/.?. 

The two electrons in the DQD are coupled by the 
Coulomb interaction, 

I e2 

where ri (r2) is the position of electron 1 (2) and n is 
the dielectric constant of the host material. In this work 
we consider the regime where the single-particle level 
spacing is the largest energy scale, in particular fvjj^ ^ 
e^/(47rKa). In this case, and assuming that the hyperfine 
and spin-orbit interactions are also weak, single-particle 
orbital excitations can be neglected. Therefore, the rel- 
evant part of the two-electron Hilbcrt space is approx- 
imately spanned by Slater determinants involving the 
Wannier orbitals $l,r- Including spin, and using second 
quantization notation where (c]j) creates an electron 
in the Wannier state ('I'i?): we define the two-electron 



basis states: 

1(2,0)5) =4^4j0), (3) 

1(0,2)5) =4^4 JO), (4) 

I(l.l)^) = ^(c[t4,i-4A)|0)' (5) 

|r+) -c[^4jO), (6) 

|ro)-^(c[,4, + 4,4,t)|0), (7) 

|r-)=c[,4,|0). (8) 



The orbital parts of the basis states with single occu- 
pancy in each well, i.e. the spin-singlet |(1, 1)5") and the 
spin-triplet |To^±), are given by 

*i(ri, r2) - [$L(ri)$R(r2) ± $R(ri)$L(r2)] , (9) 

while the orbital parts of the two states | (0,2)5) and 
1(2,0)5) with double occupation of the right and left 
wells, respectively, are given by 

*L3(ri, r2) = $L,R(ri)$L,R(r2). (10) 

The orbital functions and pj are symmetric under 
exchange of particles, and therefore must be associated 
with the antisymmetric singlet spin wave function (total 
spin 5 = 0), while ^'L is antisymmetric under exchange, 
and is associated with the symmetric triplet spin wave 
function (total spin 5=1). 

The electrostatic gates that create the potential Vdqd 
can, in addition, tune the energies of the electrons in the 
potential minima by creating an additional bias potential 
Vbias- We model this bias as a simple detuning e, which 
gives an energy difference for an electron occupying the 
left or the right dot, 

£ = ($L|eV^bias|*L) - ($R|eFbia.|$R). (11) 

In the symmetric case, e = 0, the voltages on the elec- 
trostatic gates are set so that, in the absence of electron- 
electron interactions, an electron would have the same 
energy in either well. The Coulomb repulsion, Eq.Q, 
penalizes the states |(0,2)5) and |(2,0)5) with double 
occupation of either well by an amount C/, given by 

t/=«R|t^l*L,R)- (12) 

Therefore, for a symmetric potential, e « 0, the lowest 
energy states of two electrons will be primarily comprised 
of singly occupied orbitals. When the detuning is large 
enough to overcome the on-site electron-electron repul- 
sion in one well, jej > U , the doubly occupied state with 
both electrons on the dot with lower potential becomes 
the ground state. Varying the gate voltages to increase 
the detuning from large and negative to large and pos- 
itive values then tunes the occupation numbers of the 
two dots in the ground state of the DQD through the 



sequence (2,0) — )■ (1,1) — >■ (0,2). Since the states with 
the charge configurations (2,0) and (0,2) are singlets, 
while those with the (1,1) charge configuration can be 
either singlet or triplet, the measurement of charge as a 
function of detuning can reveal the spin states. 

The strong spin-independent interactions of the elec- 
trons with the confinement potential, Vdqd, Coulomb 
repulsion C, as well as the kinetic energy are, at low en- 
ergies and in the limit of tight confinement, described 
by the matrix elements between Slater-determinant-type 
states in which the two electrons are loaded into some 
combination of the the Wannier orbitals (see Eq.(fT2|) and 
Rcf. (i2l|)): 

t = (*L3l^?,2l*R,L) - ■^iniClKn)^ (13) 

v± = (*yq*^±), (14) 

^ = (*l,r|CKl). (15) 

Here, we have used /ij* 2 — Hose + ^^DQD(r) — 14(r T ae^) 
to label the part of Hamiltonian that includes the ki- 
netic energy T and the harmonic part of the potential 
Vdqd near the dot centers. The dots are displaced by 
±a along the axis with the unit vector e^, see Fig. [T] 
Thus, /ij 2 ~ Hose describes the tunneling due to the mis- 
match between the true double-dot potential Vdqd and 
the potential of the harmonic wells [2l| . The matrix el- 
ement X describes coordinated hopping of two electrons 
from one quantum dot to the other, t is the rcnormalizcd 
single electron hopping amplitude between the two dots, 
which includes contributions of both the single-particle 
tunneling amplitude and the Coulomb interaction, and 
V+ (V_) is the Coulomb energy in the singlet (triplet) 
state with one electron in each well. 

The confinement potential Vdqd, the Coulomb inter- 
action C, and the detuning e provide the largest en- 
ergy scales in a DQD. These terms add up to the spin- 
independent Hamiltonian Hq = T + Vdqd + C + Vbias, 
which within the space of the six lowest energy dot or- 
bitals, using the basis defined in Eq. ^ - Eq. is 
represented by 

where the singlet Hamiltonian in the basis | (0,2)5"), 
1(2,0)5), 1(1,1)5) is 

/[/-£ X -V2t\ 
Hss^l X U + e -V2t , (17) 
V -V2t -V2t V+ J 

and the triplet Hamiltonian is diagonal, i?TT,o = 

In addition to the terms described above, there are 
three sources of spin-dependent interactions: Zeeman 
coupling to an external magnetic field, hyperfine coupling 
between electron spins and nuclear spins in a quantum 
dot, and the spin orbit interaction. For now we neglect 
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FIG. 1: a) double quantum dot model and the coordinate 
system. Sl and Sr denote the spin- 1/2 of the electron in the 
right and in the left quantum dot, respectively. The dots lie 
in the ^z-plane and are tunnel-coupled along the ^-direction 
(perpendicular to the z-axis). They can be detuned by the 
externally applied voltage e/e. The spin orbit field fl points 
along the z-axis, defining the first quantization axis for the 
triplet states \Tz,±) and \Tzo) (see main text). The effective 
magnetic sum-field b defines the second quantization axis for 
the triplet states |T±) and |To)- We choose the mutually 
orthogonal axes x, y, z so that b lies in the sz-plane. b) eiTect 
of detuning on the quantum dot levels. At zero detuning 
£ = 0, an electron has the same energy on the left and right 
quantum dot. For nonzero detuning, the energy of an electron 
on the left dot is e higher than the energy of an electron on 
the right dot. 



the spin orbit interaction, and analyze it in detail in the 
next section. 

The direct coupling of the electron spins to a uniform 
external magnetic field B is described by the Zeeman 
term 

i7z = -.9A*eB • (Sl + Sr) , (18) 

where g is the electron (7— factor and //g is the electron 
magnetic moment. In addition, the Fermi contact hyper- 
fine interaction between electron and nuclear spins reads 

i7„,c = 5Ih. -S,, (19) 

i 

where h^, i — L^R^ is the Overhauscr field of the quan- 
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turn dot i, given by [22 



(20) 



Here Aj is the hyperfine coupling constant for the nuclear 
species at site j, with typical size of the order of 90/ieV 
for GaAs [l^], ^'i is the electron orbital envelope wave 
function in the right (i — R) and left (i = L) dot, Rj is 
the position of the jth nucleus in the quantum dot, and 
Ij is the corresponding nuclear spin. 

Because the Zeeman and hyperfine interactions, 
Eqs. (dH]) and ([TS]). have similar forms, we combine them 
into a single effective field that acts on the electron spin 
in each dot: 
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Hz = -bL • Sl - bR • Sr 



(21) 



The effective fields bL3 = 5MbBl,r — i>-L,R include the 
contributions of the external and Overhauser fields, with 
all coupling constants absorbed in the field definitions. 
The energy levels arising from the spin-conserving Hamil- 
tonian, Eq. (|16p . along with coupling to a uniform effec- 
tive field, Eq. (PT|) with bL = b^, are shown in the left 
panel of Fig. [2] 

Below we will be interested in transitions that change 
the total spin of the pair of electrons in the DQD. To 
facilitate the discussion, we separate the total field into 
a sum- field b = (bL + bR)/2 and a difference field 6h = 
(bL - bR)/2: 



-b • (Sl + Sr) - <5b • (Sl - Sr) . (22) 



The symmetric component b conserves the magnitude 



of the total spin. 



FIG. 2: Energy levels of the double quantum dot system ob- 
tained from exact numerical diagonalization of H given in 
Eq. (|30p and plotted as a function of the detuning e mea- 
sured in units of the Coulomb on-site repulsion energy U. Of 
particular interest here are the crossings and anti-crossings 
of singlet and triplet states due to spin orbit and hyper- 
fine interactions. For Hst = (see Eq. (|42p '). i.e. van- 
ishing singlet-triplet mixing (left-hand side of plot), the pa- 
rameter values chosen are {U,t,p,X,V-,V+, b, Sby , Sh ■ e, 5h ■ 
e',ip) = (1, 0.1, 0, 0, 0.05, 0.04, 0.3, 0, 0, 0, 0). In this case, the 
singlets | (1,1)5') and | (0,2)5') anti-cross (left oval) and a fi- 
nite gap opens, whereas the singlets and triplets only cross 
(no gap). For Hst additional gaps open (right-hand 
side of plot), in particular at the lower singlet-triplet anti- 
crossing around e = e* (right oval) with an energy split- 
ting AsT that depends on magnetic field, detuning, spin 
orbit and hyperfine interactions (see main text and figures 
below). The singlet S- is a superposition of |(1,1)S) and 
1(0,2)5) (see Eq. (|37p V The parameter values chosen for 



b, (Sl + Sr; 

symmetric component 6h does not. We include the spin- 
conserving field b into the unperturbed Hamiltonian, and 
define 



0, while the anti- the right plot are ([/, t,p, X, V-,V+, b, 5by,5h ■e,Sh- e', </?) = 



TJrprp = Htt,o ~ b • (Sl + Sj 



(23) 



Below we will investigate the role of the Overhauser 
fields in causing spin transitions near a singlet-triplet 
level crossing in a two-electron DQD. While the exter- 
nal magnetic field B is a classical variable, the Over- 
hauser fields Hl.r are, in principle, quantum operators 
that involve a large number of nuclear spins, see Eq. pO]). 
Hyperfine-induced electron spin transitions may be ac- 
companied by nuclear spin flips, and the dynamical, 
quantum nature of the Overhauser field may be impor- 
tant. However, due to the large number of nuclear spins, 
the time scale for the Overhauser fields Hl.r to change 
appreciably can be much longer than the time spent near 
the avoided crossing, where spin transitions are possible. 
Thus we will treat the fields b and Sh as quasi-static 
classical variables, including a discussion of the averag- 
ing which occurs due to nuclear Larmor precession and 
statistical (thermal) fluctuations. 



(1, 0.1, 0.01, 0, 0.05, 0.04, 0.3, 0.02, 0.02, 0.01, n/2). 



III. SPIN ORBIT INTERACTION 

In addition to the external and hyperfine fields, elec- 
tron spins in a DQD are also infiuenced by orbital motion 
due to the spin-orbit interaction. Here we describe how 
the bulk spin-orbit coupling of the 2DEG is manifested in 
the confined DQD system. In GaAs quantum wells, the 
spin-orbit interaction is caused by the inversion asymme- 
try of the interface that forms the quantum well [2 4142 q 
and the inversion asymmetry of host material [27[ . With 
the 2DEG being much thinner than the lateral quantum 
dot dimensions, both spin orbit interactions are linear 
in the in-plane momenta of the confined electrons, and 
together are given by 

HsO = a {px'O-y' - Py'CTx') + P {-Px'CTx' +Py'0-v') , (24) 

where the Rashba and Dressclhaus spin orbit interaction 
constants a and /3, respectively, depend on the thickness 
and shape of the confinement in the growth direction and 
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on the material properties of the heterostructure in which 
the 2DEG is fabricated. This form of spin-orbit coupHng 
appears in a quantum well fabricated in the (001) plane 
of GaAs crystal, and the x'- and y'-axes point along the 
crystallographic directions [100] and [010], respectively. 

Within the space of low-energy single-electron orbitals 
in the DQD, the action of the spin orbit interaction can 
be expressed in terms of a spin-orbit field H, 



cLT"''c,«-h.C, 



where the field 



(25) 



(26) 



depends on the orientation of the dots with respect to 
crystallographic axes through the vector an [2^. For a 
2DEG in the (001) plane, a^j is given by 

an = (/? - a) cos Oe^^^^ + (/3 + a) sin Oe^nQ^ , (27) 

where the angle between the direction and the [110] 
crystallographic axis is denoted by 9. The matrix element 
of p^, the momentum component along the ^-direction 
that connects the two dots, see Fig. [3 is taken between 
the corresponding Wannier orbitals, and it depends on 
the envelope wave function and the double dot bind- 
ing potential. The spin orbit field accounts for the 
spin rotation when the electron hops between the dots. 
Therefore, the spin orbit interaction enables transitions 
between triplet states with single occupation of each well, 
to the singlet states with double occupation of either the 
left or the right well. 

The matrix element in Eg . (1261) c an be calculated ex- 
plicitly in a model potential (2lL 1281 ] , giving 



At I an 
3 Aso \an\ ' 



(28) 



where I is the interdot distance. The numerical prefactor 
is model-dependent, but the dependence on other param- 
eters is generic. The hopping amplitude t, and the inter- 
dot separation / depend on the geometry of the double 
dot system, whereas the spin-orbit length Ago is deter- 
mined by material properties (Rashba and Dresselhaus 
spin orbit strength) and by the orientation of the DQD 
with respect to the crystallographic axes. In particular, 
if the 2DEG lies in the (001) plane, it is given by 



1 



cos ( 



A: 



SO 



sm( 



(29) 



where X± = h/[m*{(3 ± a)] [23], with m* being the effec- 
tive band mass of the electron. In the special case /3 = 0, 
— 0, this definition reduces to the Rashba spin orbit 
length Aso |/3=e=o = ^so = h/m*a. 

One of the main goals in the following is to derive the 
dependence of the energy splitting at the anticrossing be- 
tween the lowest-energy electron spin triplet and singlet 
states (sec Fig. [21 S"- — 1+ anticrossing) in terms of this 
spin-orbit length Aso • Detailed understanding of this de- 
pendence may then be used to extract the value of Ago, 
e.g. from measurements of the singlet-triplet transition 
probability in gate-sweep experiments. 

Within the model used above for explicit calculation, 
the components of f2 are real, even in the presence 
of magnetic fields. This is due to the high symme- 
try of the ground state orbitals of the quantum dot in 
the model, and remains true even after the replacement 
PC ~^ ~ (g/c)^^ in the spin-orbit Hamiltonian H^Oi 
Eq. (P^ . However, this fact is not essential for the 
physics described below. 



IV. SINGLET-TRIPLET TRANSITIONS AND 
THE CHOICE OF SPIN QUANTIZATION AXIS 



Transitions between singlet and triplet states can be 
mediated by the spin orbit interaction or by an inhomo- 
geneous effective magnetic field (external plus hyperfine) , 
i5b. In [l^ it was shown that the transfer of angular mo- 
mentum between electrons and the nuclei strongly de- 
pends on the relative size and phase of the electron spin 
flip matrix elements induced by spin orbit interaction and 
by the difference of the Overhauscr fields in the two dots. 
Using our model of a detuned DQD, we will study these 
matrix elements in the following in detail and in partic- 
ular focus on the singlet-triplet level splitting, see Fig. [5J 
right panel. 

The homogeneous field b acts only within the 
spin triplet subspace, while the inhomogeneous field 
5h mixes singlet S = and triplet = 1 
states. Representing the total Hamiltonian in the basis 
{(1(0, 2)^), 1(2, 0)5), 1(1, l)S), 1T,+), lT,o), 1T,_))}, where 
the 2-axis is taken along rj, see Fig. [TJ we find 
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H = 



(U -e X ~V2t 

X U + e -V2t 

-V2t -^V2t V+ -V2{6b^-_i6by) 

-V2{Sb^+iSby) V-+2b^ 

iV2n i^/20 2Sb, b^V2 

V V2{6b^-iSby) 



2Sb^ V2 [Sb^ + id 
kV2 

b^V2 2b ^ 



(30) 



The Hamiltonian H , Eq. (PO)) . is the starting point for 
all of our further calculations. Our results will show the 
dependence of the singlet-triplet splitting on the param- 
eters that enter H . This Hamiltonian describes a dou- 
ble quantum dot with single orbital per dot, i.e. in the 
Hund-Mulliken approximation, and it is valid as long as 
the dot quantization energy is the largest energy scale in 
the problem. It can be applied to double quantum dots 
of various kinds, for example the gated lateral or vertical 
dots in IITV semiconductor materials, quantum dots in 
nanowircs, or self assembled quantum dots. We illustrate 
the spectrum of H in Fig. [21 for a set of parameters that 
emphasizes anticrossings of the levels. The spectrum is 
obtained by exact diagonalization of i?, and it is given 
as a function of detuning e. For other types of quantum 
dots, the parameter values would change, but the overall 
structure of the spectrum remains the same. 



V. EFFECTIVE HAMILTONIAN NEAR THE 
SINGLET-TRIPLET ANTICROSSING 

In the limit of large detuning, |e| ^ U, V+, V-, |b|, the 
ground state is a spin singlet with both electrons in ei- 
ther the left or right dot, depending on whether e > or 
e < 0. In the region of weak detuning, the ground state 
has one electron in each of the dots. For a sufficiently 
strong sum- field, |b| > U — V+, the singlet ground state 
exhibits an avoided level crossing with the lowest energy 
triplet state, i.e. the S = 1 state oriented along b, see 
Fig. [Tl at a detuning where the potential energy gained 
by the singlet's double occupancy of the lower well com- 
pensates the Zeeman energy gained by the spin-polarized 
triplet. Here, the residual splitting is determined by the 
spin non-conserving interactions. The behavior near this 
anticrossing has been the focus of many recent studies on 
the interaction of electron spins with the nuclei 
[33} . The role of spin orbit interaction has received less 
attention than that of the nuclei, and will be analyzed in 
the following sections. 

The orbital structure of the levels near the anticrossing 
is determined by the spin-independent interactions and 
by the direction and amplitude of the sum-field b. The 
singlet subspace acted on by the Hamiltonian in Eq. p7)) 
includes a state with single occupation of the two dots, 
and two states which feature double occupation of either 
the left or the right dot. Generically, the state that takes 



part in the anticrossing includes amplitudes of all three 
singlet states. However, because U ^ J, where J « 
it'^/{U-V+) ^ 0.01-0.1 mcV is the splitting between the 
lowest-energy triplet and the lowest-energy singlet state, 
and U ^ 1 meV, admixture of at least one of the singlets 
will always be suppressed at the anticrossing by a large 
energy denominator (note that t « 0.01 — 0.1 meV <C 
U-V+). 

Let us now construct the effective Hamiltonian which 
acts in the two-level subspace spanned by the levels near 
the anticrossing. First, the spin-conserving part of the 
full (6x6) Hamiltonian reads 

where the block Hss acts in the singlet subspace, Htt 
acts in the triplet subspace, and the off-diagonal blocks 
vanish due to spin conservation. 

Explicitly, the block Hss is given by Eq. pT]) in the ba- 
sis Eq. dS])-®. The triplet block Htt is given in Eq. ((23|) . 
Since U ^ t ^ X in a typical quantum dot, the state 
at the anticrossing can at most include significant con- 
tributions from two out of three basis singlets. We will 
consider the anticrossing at positive values of the detun- 
ing e (the anticrossing at negative voltage is analogous) 
so that the 1(2,0)5) state with the energy U + e is far 
detuned from the other two singlets. The remaining two 
singlets can be close in energy. In order to include the 
possibility of near degeneracy, we will introduce a mix- 
ing angle that parametrizes the hybridization of the 
1(0,2)5') and |(1,1)5) states. In the restricted subspace 
of these hybridized states, the singlet Hamiltonian reads: 

«, .£^ + .-„;2i^i_ii£ + 2*^ (32) 

where r = (Tx,Ty,Tz) is the vector of Pauli matrices. 
The pseudospin r describes the components of the an- 
ticrossing singlet state, \tz = 1) = |(1,1)5), \tz = 
— 1) — 1(0,2)5), within the approximation that we ne- 
glect the remaining |(2,0)5) component (this is valid 
when \t/U\ ^1). In this case, n = cos 2ip + sin 2ip 
is a unit vector parametrized by the mixing angle -0 that 
describes the relative size of mixing and splitting of 
eigenstates. Note that the y-component of n vanishes due 
to the choice of phases in the quantum dot ground states, 
which guarantees that the spin-independent hopping ma- 
trix element t is purely real. We remark again that in our 
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DQD set-up the hopping matrix element t stays real even 



in the presence of magnetic fields, j21| . The mixing angle 
of the doubly and singly occupied states at the S = an- 
ticrossing, |(0,2)S') and Kljl)^), respectively, is defined 

by 



cos 2?/' = 
sin 2'ip = 



U-V+-e 



{U-V+- ef + 8^2 
^{U~V+~ ef + 8^2 



(33) 
(34) 



In the basis of eigenstates of Eq. (|32)) the singlet Hamil- 
tonian _ffss is given by 



Hi 



ss 



U 

X cos ij) ~ 
—X sm^j - 



\/2t sin i\i 
- \/2t cos Tp 



X cos i/j — sin -0 — X svnij} — y/2t cos ijj 
Es+ 



(35) 







Es- 



where 



E. 



U -e + V+ 



s± 



F+)2/4 + 2t2 (36) 



are the eigenvalues of iJgs in the spin-conserving sector. 
The basis vectors used here are the far-detuned singlet 
1(2,0)5) (Eq. (P) and the singlets \S±) defined by 



l^+) 
\S-) 



sin?M(l,l)S') - 
cosVKl,!)^) 



cos V'l (0,2)5), 
- sin V'l (0,2)5). 



(37) 
(38) 



In the limit U ^ t, \S±) become eigenstates of iJss with 
energies Es± given in Eq. (|36|) . 

While the mixing of orbital states belonging to sin- 
glets does not affect the triplet Hamiltonian Htt, it 
will change the form of coupling between the singlet and 
triplet states near the anticrossing. In the following, we 
first diagonalize the triplet sector in order to find the ex- 
plicit form of the triplet state at the anticrossing, and 
then find the effective Hamiltonian of the singlet-triplet 
coupling. 

The triplets are Zeeman split by the sum-field b. We 
have chosen the z-axis of spin quantization so that the 
spin-orbit interaction couples |(0,2)5) and |(2,0)5) to 
the \S = l,Sz = 1) state. We will now diagonalize the 
triplet part of the spin-conserving Hamiltonian, given by 



H' 



TT 



V-+2b 



■ COS (y5 

^sin^ 




^sm^ 


-ijsin^ 





sin ip 



1 

cos ifl 



(39) 



where we have used h = |b|, cos(p = 6z/6, and sin(p = 
hxl^ (see Fig.[T]). The unitary transformation HttU^ 
that diagonalizes Htt is 

/ 



cos 
1 



2 £_ 
2 

sin tp 

2 V 




(40) 



We denote the basis states by |T+), \To), and |T_), where, 
now, the quantization axis is given by the sum-field b. 

With the diagonalization of the triplet block of the 
spin-conserving Hamiltonian, and the preceding ap- 
proximate diagonalization of the spin-conserving singlet 
Hamiltonian, Eq. psp . we are able to describe the spin- 
conserving interaction near the anticrossing in a conve- 
nient form. We will use |5_) and |T+) as basis vectors, 
and study the effective Hamiltonian in the vicinity of the 
anticrossing that can emerge from spin non-conserving 
interactions. 

The full Hamiltonian is 



H — Hsc + Hso 
and reads in block form 



Sh ■ (Si 



H 



(41) 



(42) 



The diagonal blocks Hss and Htt give the spin- 
conserving part, denoted by Hsc, while the off-diagonal 
blocks HsT and Hts = -f^sT induce singlet-triplet transi- 
tions. In the basis {1(2,0)8) ,\S+) ,\S-) ,\T+) ,\To) ,\T_)), 
the diagonal blocks take a simple form. The singlet block 
is _ffss given in Eq. The triplet block Htt is diag- 

onal and reads 



Htt = diag (K_ - |b|, V.,V. + |b|) , 



(43) 



where we used that |r+) is the lowest-energy triplet. 

The effective Hamiltonian near the anticrossing is de- 
termined by the spin-conserving terms t, U, X, e, b, and 
by the spin non-conserving terms 17, arising from spin- 
orbit coupling, and 6h, the effective difference field. All 
these interactions can be treated perturbatively in dots 
with weak overlap of the orbitals. The off-diagonal terms 

are denoted by Hst = i^fx ' u wt 

where 



Hg^, and Hts — -^sT' 
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ttSO 
-"ST 



— sin(y3 -~^/2 cos ip siiKf 

ifl I cos "4! sin If \/2 cos "0 cos (ys — cos ip sin ip 
— sin ■(/) sin — \/2 sin t/; cos (/3 sin?/'sin(y5 



(44) 



5b 

ST 





\/2 {iSby + 5h ■ e') sin ip 2(^b-e)sinV' \/2 {i6by — Sh ■ e') simp 



\/2 {i6by + 6h ■ e') cos ip 2((5b-e)cosV' V2{iSby 



(45) 



5h 



I cos 



where the unit vector 

e = ex sin ip + ez cos 



(46) 



points in the direction of the homogeneous field b, and 
the vector 



-e^ cos (p + Bz sm ip 



(47) 



hes in the xz-plane, which contains Jl and b, and points 
in the direction normal to b. 

In the vicinity of the anticrossing, the DQD behaves 
as an effective two-level system, with the dynamics de- 
scribed by an effective Hamiltonian denoted by Her- Up 
to first order in spin non-conserving interactions and, af- 
ter neglecting the high-energy state |(2, 0)5"), we find 



(1) _ j Es- H34 
Hi3 Et+ 



WP> = 



(48) 



where H34 = {S-\H\T^) is the matrix element of H, 
Eq. (|42t . between the anticrossing states, and Et+ = 
V--b. 



VI. SINGLET-TRIPLET SPLITTING AT THE 
ANTICROSSING 

The singlet-triplet splitting Ast at the S- — T+ anti- 
crossing, see Fig. [51 can be accessed in the spin blockade 
regime by transport measurements or by charge sensing. 
This splitting gives valuable information about the prop- 
erties of the quantum dots and the nuclear polarization. 
We shall derive now explicit expressions for Agx in terms 
of the experimentally relevant quantities e, b, (5b, and Jl. 
We proceed by perturbation expansion in -ffsT, t/U, and 
X/U . First, we focus on the first-order contributions and 
afterward address the higher order corrections which be- 
come relevant around the points where the leading contri- 
butions can be tuned to zero by the control parameters. 

As the DQD detuning e is varied with other parameters 
held fixed, there is a special value e* where the energy 
Es- of the lowest singlet, Eq. (p6| . becomes equal to the 
energy Et+ ^V^-boi the triplet \T+), Eq. The 
detuning at which this crossing occurs is controlled by 
the amplitude b of the sum-field, as well as the tunnel 



coupling t. For the unperturbed case, described by iJo, 
e* is the solution of the equation 



Es-{e*) = V- 



(49) 



When the spin non-conserving interactions are taken into 
account, the crossing of singlet [S — 0) and triplet {S = 
1) states is avoided due to state mixing (hybridization). 

Up to first order in -ffsT, t/U , and X/U , the splitting 
follows from Eq. ([^S]) and reads 



Ast(£, ^ </?) = 2, /|i/34|2 + (Ss- - ^T+)2 



(50) 



For a fixed value of 6, the splitting attains its minimum 
value Agrp at e*, 

A^T = AsT(e*(^),6,<^) = min Ast (£,&,'/'), (51) 

e 

where e* implicitly depends on 6, as well as other DQD 
parameters. From Eqs. (|48)) . (|44)) . and (|45|) . the splitting 
is equal to 2|i734|, 



A, 



ST 



-iVt sin p sin tp + \pi ((5b • e' + ibby) cos ip 



(52) 

Note that Agrp contains contributions from both the spin- 
orbit coupling and the difference field (5b. The relative 
importance of each of the two terms depends on the de- 
tuning e* through the mixing angle ip, as well as on the 
geometry through the angle ip between the effective field 
b and the spin-orbit field fl. When varying the detuning 
e from large to small values, tp decreases from -ip ^ it/2 
at strong detuning, £ [/ — V+, to 1/; « at |£| <C t. For 
a mixing angle tp sa 7r/2, the contribution to Agrp coming 
from the spin orbit interaction dominates the one from 
the difference-field, and vice versa for ■)/' ~ 0. 

Reaching the regime requires weak magnetic 

fields, i.e. b -^t, and in this case the energy splittings be- 
tween the triplet states are not large enough to make use 
of the simple model for a two-level anticrossing. On the 
other hand, reaching the regime with -0 sa 7r/2 requires 
that the detuning e* at which \S-) and |r+) anticross 
is far away from £12, the detuning at the anticrossing of 
1(1,1)5') and 1(0,2)5') singlets, cf. Fig. [2] The width of 
the 1(1,1)5) — 1(0,2)5) anticrossing is of the order t, so 
the requirement is |£* — £12] ^ t. Therefore, the Zee- 
man energy of |T-|_) must be larger than i, so q^ibB ^ i. 
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which gives B ^ 0.2 T for typical values t 10 /icV and 
l5l=0.4. 

These considerations show that, at least in principle, 
the relative strengths of the spin-orbit and hyperfine con- 
tributions to the singlet-triplet coupling can be tuned 
through a wide range of values using a combination of 
gate voltages and magnetic field strength and direction. 
What situation do we expect for typical GaAs dots? Us- 
ing a value of 3-5 mT for the random hyperfine field (see 
e.g. Ref. m and references therein), and an electron g- 
factor l^l = 0.4, we estimate \6b\ « 70-120 neV. For the 
spin-orbit coupling il, sec Eq. ipS)) . using t ~ 10 fieV, an 
interdot separation I = 50 nm, and a spin-orbit length 
Aso in the range 6-30 /um (see e.g. Refs.[l0ll33|). we find 
|0| « 20 — 110 neV. Parameters may vary from device 
to device, but it appears that the spin-orbit and hyper- 
fine couplings are generally of similar orders of magni- 
tude, with the spin-orbit coupling typically a few times 
weaker. Thus adjustments of the matrix elements over a 
reasonable range of ip may be sufficient to explore both 
the hyperfine and spin-orbit dominated regimes. Similar 
analysis can be performed for devices in other materi- 
als, such as In As or InSb nanowires, where the natural 
balance between hyperfine and spin-orbit couplings may 
shift. 




FIG. 3: Singlet-triplet level splitting Ag^ = Ast {e* ,b,ifi), 
Eq. (|52|l . for 5h — (no nuclear field), as a function of the 
detuning e* at which the crossing occurs (cf. Fig. |2]) and of 
the angle (p between the spin-orbit field fl and the magnetic 
sum-field b (cf. Fig. [l]). Parameters used for this plot are 
U = l,t = 0.01, V+ = 0.75, (5b = 0. 



A. Singlet-triplet splitting Ag^ for 5b = 

Let us first consider the special case of vanishing dif- 
ference field, Sh ~ 0, and finite uniform field, b 7^ 0. In 
this case, the splitting depends not only on ft but also on 
the detuning at the minimal splitting, e*, which itself is 
implicitly determined by the sum-field amplitude b. The 
geometry of the system enters through the angle ip be- 
tween b and the spin-orbit field $7. In Fig. [3] we show a 
plot of the level splitting Ast {£* ,b,ip), given in Eq. ((5^ 
as a function of its variables and for 6h = 0. 

At any fixed angle f ^ 0, Agrp shows a dependence on 
the detuning e* at the anticrossing due to the mixing of 
1(0,2)5'), which is coupled to |T+) via the spin-orbit in- 
teraction, and 1(1, 1)S), which does not couple to \T^) via 
spin-orbit coupling in the first order, see FigU) At large 
values of detuning e* U — V+, the sphtting reaches 
a saturation value 2Qsm(p. For typical GaAs quantum 
dots, reaching this regime requires strong magnetic fields 
of B ^ t. At lower values of the detuning e < J7 — V+, 
the mixing of the singlet 1(1,1)5') becomes significant, 
and the spin-orbit coupling value f2 cannot be read off 
directly from the splitting. 

The maximal splitting Agrp caused by spin-orbit inter- 
action is Agrp = 20|sin(^|, for ip = 7r/2. From Eq. (pS]) 
and Eq. we find that is set by the material prop- 
erties (Dresselhaus (/3) and Rashba (a) spin orbit inter- 
actions) and the geometry of the dots. Assuming that 
the magnetic field is strong enough to separate the |To) 
and |r_) states from the anticrossing, U ^ b ^ t ^ |5b|. 



the maximal splitting is (| sin -01 = 1) 
At I 

Al^ = - — \,m^\, {5h = 0), (53) 

where I is the interdot distance. The numeric factor (of 
order unity) is non-universal and depends on the spe- 
cific dot geometry. Formula (|53p is one of the main 
results of this paper. It provides a simple but useful 
relation between quantities that can be determined ex- 
perimentally, such as Agrp, t, /, and ip, and a quantity 
of interest - the spin orbit length Aso- This relation 
could allow the strength of spin-orbit coupling to be 
measured experimentally [l3, Hj], though the geometry 
and detuning-depcndence must be carefully taken into 
account in order to obtain an accurate estimate. 

Let us remark here briefiy on the special case of zero 
detuning, i.e. e = 0, and weak magnetic fields. In this 
case, the splitting is not described by our calculations 
which require sufficiently large separation of the triplets 
in energy. However, in a slightly different system ~ a 
single quantum dot containing two electrons - singlet- 
triplet coupling which is forbidden by time-reversal sym- 
metry can be generated by applying a magnetic field, 
Agrp w (aB/Aso)-E'z [l^- This case cannot be recovered 
from our DQD model with one orbital per site. Indeed, 
here we have seen that in weak fields, b <^ t, the cou- 
pling of the two states with single occupation in each 
well 1(1,1)5) and |r+) due to the spin-orbit interaction 
involves doubly occupied states which are higher in en- 
ergy due to the on-site repulsion. On the other hand, for 
a pair of electrons in a single quantum dot, the on-site 
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£/U 

FIG. 4: Singlet-triplet level splitting Agx = ^ST{e* ,b,ip), 
Eq. (|52|l . for (5b = (no nuclear polarization), as a function 
of the detuning e* at which the crossing occurs (cf. Fig. [2]) 
and of the angle ifi between the spin-orbit field fl and the 
magnetic sum-field b (cf. Fig. [T]). At small detuning e* < t, 
the splitting becomes rather small, while it saturates at large 
detuning e* > U. The saturation value is oc |sinyi|, shown 
for (/3 = 7r/4 (full line), if = Sn/8 (dashed line) and ip — 5n/8 
(dashed-dotted line). Parameters used for this plot are U = 1, 
t = 0.01, V+ = 0.75, 5h = 0. 



repulsion is approximately the same for both states, sin- 
glet and triplet. We note that for DQDs in weak fields, 
the Zeeman energy Ez, occurring in the splitting for a 
single dot , gets replaced by the exchange energy J if 
the mixing of triplets due to (5b 7^ is neglected. 




FIG. 5: The same plot of the singlet-triplet splitting as 
in Fig. [3] except for finite nuclear polarization chosen to 
be Sh = (-0.0006,0.0008,0.0012). In the strong detuning 
regime, e* > U — V+ on the right-hand side of the plot, the 
splitting is determined by spin orbit interaction, uso > '^hf 
and resembles the same area in Fig. [3] In the weak detuning 
regime, e* < U — V+ on the left-hand side of the plot, the 
hyperfine interaction increases the splitting. The regime with 
similar strengths of the interactions, vso ^'hf can be iden- 
tified in the region e* ~ [/ — V+ = 0.25?7. For e* on the left- 
hand side of the vertical dashed line, and ip ~ 0,tt, the split- 
ting is dominated by hyperfine interaction vhf > vso- Near 
the dotted line, and for the angles tp ~ n/2, 3tt/2 the sizes of 
spin orbit and hyperfine interactions are similar vso ^ vhf- 



B. Singlet-triplet splitting Ag^ for (5b 7^ 

In addition to spin-orbit coupling of the anticrossing 
triplet to 1(0,2)5), the anticrossing triplet is coupled to 
the singlet |(1, 1)5) through the difference- field 6h. The 
previous considerations show that the contributions from 
the difference-field Sh to the splitting cannot be neglected 
for angles (p ^ 0, tt, or for field strengths where e* < 
[/ — V+, which is often be the case. Therefore, we now 
discuss the splitting in the presence of both, the spin- 
orbit field ^l and 5h. The splitting Agrp as a function of 
detuning e* and direction of b is shown in Fig. [5] 

When both sources of splitting are present, generically, 
the gap Agrp remains open. The spin-orbit contribu- 
tion to Eq. (|52p is always purely imaginary, while the 
(5b-contribution has both a real part, coming from the 
component lying in the a;2;-planc, and an imaginary part, 
coming from the perpendicular component Sby. For any 
fixed value of the spin-orbit coupling strength, closing the 
gap would require fine tuning of i5b, both in amplitude 
and direction. As a function of the direction of b the con- 
tributions to Agrp compete, and, in addition, the relative 
size of the competing terms will change as a function 



of s* . Indeed, the spin-orbit term is strongest at large 
detuning e* ^ U — V+, while the difference-field term 
becomes significant at small detuning s* ^ t. This com- 
petition affects the form of the (/j-dependent splitting, see 
FigEl 

In the limit |f2sin'(/'| ^ | (5b cos ■0 1 , the splitting Agrp 
is caused mostly by the inhomogeneous field. In this 
case, the splitting is proportional to the size of compo- 
nent (5b = (5b — e((5b • e) of (5b which is normal to the 
homogeneous field b = |b|e. With the leading spin-orbit 
coupling correction, the splitting is [see Eq. ([5^ ] 



^ST 



2\/2|(5b_L cosV'l 



2^5hy sin il) cos ip 
|(5bx cos '01 



(54) 



C. Measuring the singlet-triplet coupling 

The singlet-triplet coupling Agrp is manifested experi- 
mentally, for example, in the spin flip probability when 
the system is taken through the level crossing during a 
time-dependent gate sweep (35| . In such experiments, 
the system is initialized to its ground state at large e, the 
(0, 2) singlet. When e is then ramped to take the sys- 
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FIG. 6: First order singlet-triplet splitting Ag^ — 
AsT(e*, 6, V^), Eq. (|52|) . for finite nuclear polarization 5b 7^ 0, 
plotted as function of the angle ip (cf. Fig. [TJ, for various 
detunings e* (cf. Fig. [5} . The parameters used for the plot 
are [/ = 1, t = 0.01, V+ = 0.75, V- = 0.74, J7 = 0.005, 
and Sh = (-0.0006,0.0008,0.0012). The curves correspond 
to e* = 0.1 (full line), e* = 0.2 (dashed), e* = 0.4 (dashed- 
dotted). In the strong detuning regime (dashed-dotted line) 
the angular dependence reflects the | sin dependence of the 
spin-orbit term (see Eq. (|53|) . The regime of weaker detuning 
(full and dashed line) shows the hyperfine effects. 




0.4 0.6 
(p/27l 

FIG. 7: Square of the splitting (Ast), Eq. ([52|l . averaged 
over Gaussian fluctuations of 5h with zero mean and standard 
deviation a. The plots show the dependence of (Ast) on 
the angle ip for various mixing angles tp. We have assumed 
isotropic Gaussian fluctuations with a standard deviation a. 
Plots are for the values a/Q — 0.1, 0.5, 1, 2 and illustrate the 
effects of various strengths of fluctuations. The curves are 
found by numerical averaging over the fields 5b. 



tern through the singlet-triplet crossing, the two-electron 
spin state may change, with a probability determined by 
a combination of the coupling Ag^p and the sweep rate. 
The final spin state can then be read out by quickly ramp- 
ing back to large e, where the singlet and triplet states 
have discernibly different charge distributions, which can 
be detected by a nearby charge sensor. 

Even with single-shot spin detection determin- 
ing the spin flip probability requires building up statis- 
tics over many experimental runs. Within each run, the 
parameters in Eq. (j52p may be considered fixed. How- 
ever, the hyperfine field components are in general a pri- 
ori unknown: under typical experimental conditions, the 
temperature is high compared with all intrinsic energy 
scales within the nuclear spin system, and the equilib- 
rium state is nearly completely random. Depending on 
the measurement timescale, the nuclear fields on subse- 
quent experimental runs may either remain constant or 
may change. While the correlation time for the longi- 
tudinal component of the nuclear field (parallel to the 
external field) may be quite long, the transverse compo- 
nents change on the timescale of nuclear Larmor preces- 
sion, which for moderate fields of a few hundred millitesla 
can reach the sub-microsecond timescale. The coherence 
time associated with this precession may reach several 
hundred microseconds to one millisecond. 

Let P(Agrp) be the probability that the system makes a 
transition to the triplet state in a single sweep, when the 
value of Agrp is specified. In an experiment where mea- 
surements of the singlet and triplet fractions are averaged 



over a time long compared to all nuclear spin relaxation 
times, one obtains an averaged probability (P(Agrp)), 
where (A) denotes the mean value of quantity A, av- 
eraging over a Gaussian distribution of 6h, while other 
parameters such as B, t, 6, (j) and the sweep rate are held 
fixed. If the measurements are averaged over a shorter 
period, which is long compared to the time for phase re- 
laxation of the nuclear spins, but short compared to the 
longitudinal relaxation times, then the Gaussian average 
should be taken only over the transverse components of 
6h, while the component parallel to the applied magnetic 
field is held fixed. 

When the sweep rate through the S-T transition is 
rapid, the probability P(Agrp) should be proportional to 
(Agrp)^, so an average value of P(Agrp) will measure the 
mean value of (Ag^)^. For lower values of the sweep 
rate, P will have corrections due to (Ag^p)^, etc.. There- 
fore, measurements of the averaged value (P(Agrp)) for a 
wide range of sweep rates should, in principle, yield av- 
erage values of all powers of (Agrj,)^, and thus allow one 
to deduce the probability distribution for (Agrp)^. Here 
we concentrate on the mean value of (Agrp)^ , and dis- 
cuss predictions for this mean value as a function of the 
parameters B, (j), and t. 

We illustrate the dependence of ((Agrp)'^) on the angle 
if and mixing ?/' in Fig. [7] and Fig. [51 The dependence of 
the splitting on the angle f, inherent in the non-averaged 
splitting, see Eq. (j52p . remains visible when the splitting 
is dominated by spin-orbit interaction. As expected, the 
dependence of the splitting on the angle (p, Fig. [71 is most 
visible in the case of weak fluctuations of Sh, i.e. for weak 
hyperfine coupling. In addition, the angular dependence 
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FIG. 8: Average square of the splitting (Agx), in the fluctu- 
ating nuclear field 5b. The parameters are chosen as in Fig.[7l 
and we illustrate the dependence on the mixing angle tp. 



is more pronounced for larger mixing angles, since the 
spin-orbit induced splitting depends on flsmip, Fig. [5] 

The fluctuating difference field, besides changing the 
average ((Agrp)^), introduces noise in the splitting. We 
find that the standard deviation cta of the spfitting, in 
the limit of weak fluctuations irJsint/ssin-i/)! ^ cr is 

o-A.G = 2o-| coaip\, (55) 

so that it also shows dependence on the mixing angle ip. 

Fluctuations in Sh smear the splitting at the anticross- 
ing Agrp. The average value and noise in the splitting 
can be used to measure the strengths of spin-orbit cou- 
pling ri and the hyperfine field 5h. The relative size of 
fluctuations in Agrp, function of ip, has minima at 
ip « 7r/2,37r/2. 

D. Higher order corrections 

So far, our analysis of the \S-) - |T+) anticrossing was 
based on the assumption that the largest contribution to 



the splitting Agrp, Eq. (|52|) . results from the direct cou- 
pling of the two states via Hg^ and H^^. However, if the 
detuning e* is not large enough to make the influence of 
the levels that are energetically further away from the an- 
ticrossing completely negligible, higher order terms that 
describe virtual transitions to such higher levels and thus 
involve more than one transition between the singlets and 
the triplets become important. 

To study this regime, we derive an effective Hamil- 
tonian in the vicinity of the anticrossing by a second 
order Schrieffer- Wolff (SW) transformation [sl, [13 ■ We 
divide the Hilbert space of the DQD into a relevant 
part which includes the anticrossing states, and an aux- 
iliary part which contains the remaining 4 states. The 
time-independent perturbation series is then performed 
in powers of the spin-non-conserving interactions. The 
spin-conserving Hamiltonian Hgc, Eq. (|3ip . is taken as 
the unperturbed part, while H — Hgc is the perturbation. 
In reordering the basis, we choose the crossing states jS*-) 
and |T+) of Hgc to be the first two basis states. Then, 
the Hamiltonian has a block-diagonal form denoted by 



where A is a 2 x 2 matrix that describes the an- 
ticrossing states, i3 is a 4 X 4 matrix of the states 
with energies far from the anticrossing and the 2 x 
4 matrix C represents the coupling between of the 
subspaces controlled by A and B. In the basis 
(|S'_),|r+),|ro),|r_), 1(2,0)5), |5+)), with the two anti- 
crossing level at the positions 1 and 2, we can read off 
the blocks from Eq. ((55)) : 



A 



Es- 

iilg sin tp + ^/2S- cos ip 



-ifls sin-0 -I- \/2S+ cosi/j 



E' 



(57) 



B 



( V- isfklc. -iV2ficCosV' + 2(5sin'i/> \ 

V- +h — ifis ifis cos t/i — a/2i5+ sin 

~i\f%lc i^s U + e ~ \/2t sin ip — X cos ip 

\ i^/2flcCos^ + 26 s'm^p —iflg cos ^ — ^/2S^ sin ip — \/2t sin ip — X cos ip Es+ ) 



, (58) 



C ^ 



i-\/2fic sin?/; -I- 2(5cos'0 if2s sin?/; — \/2<5- cost/' ~ \f2t cos ip + Xsni'xp 



-Klg cos -ip sin -0 



(59) 
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where we have used the abbreviations S = Sh ■ e, 5± = 
5h ■ e' ± iSby, fls = ^sinip, Vl^ = Qcosip, and the unit 
vectors e' and e are defined in Eq. (|47| and Eq. (|46| . 
respectively. 

As a result of the SW transformation on Eq. (|56p , the 
off-diagonal block C is eliminated up to second order in C 

and the transformed block A A + SA becomes then 
the Hamiltonian of an effective two-level system. There- 



fore, the first-order Hamiltonian Her from Eq. (|48|) be- 
comes modified by second order terms, Her Her +5 A, 
where 6 A is the second order correction to A. The diag- 
onal matrix elements SAn and 6A22 describe the renor- 
malization of the energy levels, and their effect is to shift 
the detuning e* at which the anticrossing occurs. The 
explicit expressions of these corrections are 



SAu 

SA22 
M12 



Es- - V- 



2(17cSin?/') +4r cos^i' 



1 



£^5- -V--b 



Es--U-e 
2 ((5b • e')^ cos^ 7/1 + {^/25by cos-0 + sin^i/; 
1 



\f2t cos V' - X sin V-') + 



{^s cos i/) -I- \f28hy sin t/i^ -1-2 [Sh • e')^ sin^ 



-VL. 



(fit cos V' + X sin t/i 



1 



1 



Es- - U -e Et+ -U -e 
I 



(60) 

(61) 
(62) 



In experiments that probe the electron spin dynam- 2\H34^ + SAi2\. Thus, up to second order, the splitting at 
ics, the most important terms are the off-diagonal ones, the anticrossing becomes 
SA12 = SA21. They lead to a modification of the first 

order singlet-triplet splitting Eq. ([STj). i.e. Agrp 



Agrp — 2 



—iQ sinip 



sin?/' 



—= cos ip sm ip 

v2 2 



Es_ - U -e* Et+ - U -£* 



V2 ((5b • e' -t- i8hy) cos ip ■ 

m 



We see now that the new correction terms in Agrp become 
significant for weak detuning, when ip < tt/2, because 
the spin-independent tunneling contribution, which is 
oc cosip, can alter the first-order result, which is oc simp. 
We compare the splitting in the second order, Eq. (|63l). 
with the first order splitting and the result of exact nu- 
merical diagonalization of H, Eq. (pO|) . in Fig. [9l The 
higher order corrections are small, but they do become 
significant for the magnetic field normal to the spin or- 
bit parameter, (p = 7r/2,37r/2, due to stronger effective 
strength sin (p of spin-orbit coupling. For other consid- 
ered values of detuning, e* = O.IU and e* = 0.2U, the 
change of splitting is smaller than for the e* = OAU case. 

The limit ip — 7r/2 requires strong magnetic fields, 
B '--^ 1 T for a typical GaAs DQD. It is reasonable to as- 
sume that the experiments can be performed both in this 
limit and away from it, so that the dependence of Agrp 
on ip can be probed. In materials with larger g-factors 
such as InAs, InSb, SiGe, the limit is reached at lower 
fields. In addition, we have obtained similar results for a 
model DQD with t = OAU, fl = OAt, \Sh\ w 0.5^, that 



describes a smaller DQD with more pronounced hopping. 



VII. RATIO OF SPIN-ORBIT AND HYPERFINE 
TERMS 

As pointed out before, the spin non-conserving Hamil- 
tonian in the vicinity of the anticrossing can be used to 
get experimental access to the spin orbit interaction and 
nuclear polarization in the difference-field (5b. Being a 
function of controllable parameters e and b, this Hamilto- 
nian can be altered by applying voltages to electrodes in 
the vicinity of the quantum dots, adjusting the strength 
of an external magnetic field, or changing the direction 
of the field. 

The effect of the competition between spin orbit and 
hyperfine induced spin flips on the efficiency of angular 
momentum transfer between electron and nuclear spins 
was recently studied theoretically, both in the context 
of dc transport experiments [isj and in the context of 
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— first order 

— second order 
— numeric 



e = 0.4 U 



0.2 



0.4 0.6 
(p/ 2 TT 



0.8 



FIG. 9: Comparison of the splitting Ag^^ obtained from the 
perturbation in first order (full line), Eq. (|52|l . and in the 
second order (dashed line), Eq. (|63|) to the exact numerical 
result (dashed-dotted line), obtained by the direct numerical 
diagonalization of H , Eq. (|3Up . The plots show the splitting 
Agx as a function of the angle ifi for the anticrossing at the 
detuning e* = 0.4C7. The parameters used in this plot are 
U = 1, t = 0.02, V+ = 0.75, V- = 0.74, Q, = 0.005, and 
5h = (-0.0006, 0.0008, 0.0012). 



gate-sweep experiments [l^ [1^ . References [l^ and [l^ 
revealed striking sensitivities of the polarization transfer 
efficiency on the ratio of spin orbit and hyperfine cou- 
pling strengths. In those works, the coupling strengths 
were treated as phenomenological parameters. Here we 
provide explicit expressions for them, and discuss how 
they can be tuned. 

Following the notation of Ref. [l^ we write 



= vso + e*''wHF, 



(64) 



where wgo and whf stand for the transitions caused by 
spin orbit and hyperfine interactions, respectively. Com- 
paring with Eq. (|48|) . we can identify with Then, 
adjusting the overall phase to make the spin-orbit part 
Vso real, we identify in lowest order 



Vso ~ I sin sin -0 1 



WHF = I cosV'l Y (<5b • e') +(56,2, 
"-(5b e'" 



•d = arctan 



(5fe„ 



(65) 
(66) 

(67) 



The explicit expression for d shows that the phase of the 
matrix element can be adjusted not only by changing 
the direction of (5b. but also by rotating the external 
magnetic field, which changes e', and also controls the 
effective spin-orbit coupling strength. 

Our results show that, in principle, it is possible to 
switch between the two regimes of hyperfine-dominated 
and spin-orbit-dominatcd behavior, by changing the ex- 
ternal magnetic field strength and direction. In the 
strong field regime, with sum-field b being large, e* is also 



large, and thus the spin-orbit terms become dominant (V' 
approaches 7r/2). This behavior is illustrated in Fig. [S] 
Note that the large values of e* require strong h fields. 
On the other hand, switching to the regime whf > wso is 
always possible by rotating the direction of the magnetic 
field so that it coincides with ±Jl/|r2| giving (y9 w 0. In 
this case, the term vso is negligible and uhf dominates 
the splitting. Higher-order corrections to the effective 
Hamiltonian at the anticrossing point do not alter this 
basic picture of the splitting, but they do change the val- 
ues of the parameters ■0 and Lp at which the switching 
occurs. 

The switching between the regimes dominated either 
by spin orbit or by hyperfine interactions can potentially 
be achieved as follows. For the regime vso > 'I'hf, the 
sum-field b should point along the spin orbit field il, see 
Eqs. (|26p and (|27|) . in order to maximize |sin(^|. Also, 
the applied field should be as strong as possible, in order 
to maximize the amplitude of the singlet |(0, 2)5*) (con- 
tributing to the anticrossing singlet \S-), see Eq. ((57|) ). 
On the other hand, the opposite regime, whf > vsoi can 
be reached by orienting b along $7, and thus reaching 
sin</7 = 0. If sin(/5 = cannot be achieved, i;so can be 
reduced by decreasing b and thereby increasing the am- 
plitude of the singlet state |(1, 1)5*) in the |5'_)-singlet at 
the anticrossing. 



VIII. CONCLUSIONS 

We have derived an effective two level Hamiltonian iJcr 
for a detuned two-electron double quantum dot in an 
external magnetic field. Our effective Hamiltonian de- 
scribes the dynamics of the electron spins for the values 
of detuning e ~ e* close to the anticrossing of the lowest 
energy 5 = and the lowest energy S = \ state. We have 
shown how -ffcr can be used in the interpretation of ex- 
periments that probe electron spin interactions by charge 
sensing and transport in the Coulomb blockade regime. 
The dependence of Il„ on the detuning and magnetic 
fields can also be used to switch the spin dynamics in a 
double quantum dot between the spin-orbit dominated 
regime, and the hyperfine-dominated regime. 

The spin dynamics at the anticrossing is governed by 
the spin-orbit and nuclear hyperfine interactions. In a 
double quantum dot, these two interactions act differ- 
ently on the orbital electronic states. On one hand, the 
spin-orbit interaction causes hopping of an electron be- 
tween the quantum dots accompanied by a spin rotation, 
thus changing the occupation of the quantum dots. On 
the other hand, the nuclear hyperfine interaction acts as 
an inhomogeneous magnetic field, and causes spin rota- 
tions that are local to the dots, leaving the charge state 
unchanged. Due to this distinction, the detuning e con- 
trols the relative strength of the two interactions in i/cr, 
in addition to the ratio |r2|/|5b|, or |il|/(T. In the limit 
of detuning much stronger than the on-site repulsion of 
the dots, e ^ [/. i/cr describes mostly the spin-orbit in- 
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teraction, with negligible hyperfine effects. In the case 
of weaker detuning, the effective hyperfine interactions 
can be of the size comparable to the effective spin orbit 
interactions. 

In addition, we find that the orientations of both the 
double quantum dot and the external magnetic field, de- 
scribed in the Hcv by the spin-orbit field and the sum 
field b, infiuence the effective spin orbit interaction. In 
particular, by having b pointing along VI, we can sup- 
press the spin orbit effects completely (in leading order). 

The splitting of the anticrossing states is accessible 
to experiments. It can be calculated from ifcr, and we 
find the dependence of this splitting on detuning and the 
strength and direction of the sum field, AsT(e, (p). Of 
particular interest is the splitting of levels at the anti- 
crossing. We calculate this quantity, Agrp(e*, (/j), as a 
function of the detuning at the anticrossing point, £*, 
and the orientation of the sum field, given by the angle 
Lp. Both the spin orbit interaction strength and the in- 
homogeneity in the hyperfine coupling can be deduced 
by measuring the splitting and using our formulas for 
A*T(e*,^). 



The relative strength of the spin orbit and hyper- 
fine terms in Her has a profound effect on the coupled 
dynamics of electron and nuclear spins. The value of 
the average angular momentum transfer to nuclear spins 
as an electron tunnels through a spin-blockaded DQD 
changes sharply as the interaction goes from the spin- 
orbit-dominated to the hyperfine-dominatcd regime. The 
spin orbit interaction is dominant in the limit of strong 
detuning e* 3> U — V+. The regime dominated by nu- 
clear hyperfine interaction is reached when the detuning 
is weaker e* t and the orientation of the sum field 
is along Vl. Using the dependencies of the matrix ele- 
ments on gate voltages and magnetic field strength and 
orientation, it may be possible to tune between these two 
regimes in situ, thus enabling experiments to study their 
sensitive competition. 
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